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ABSTRACT  'j  ^ 

One  of  the  more  M.mportant  phenomena  associated  with  the  upward 
motion  of  a  gas  bubble  from  an  underwater  explosion  is  the  significant 
departure  from  sphericity  near  the  times  of  the  minimum  bubble  radius* 
Neglecting  this  change  in  shape  results  in  the  prediction  of  a  much 
faster  upward  velocity  than  actually  occurs.  -The  inclusion  of  this 
effect  in  the  equations  of  motion  has  been  exceedingly  difficult  be¬ 
cause  of  the  large  magnitude  of  the  departure  from  sphericity.  '  .  /  / 

In  this  work,  the  shape  of  the  bubble-ls  described'  by  an  ellipsoid 
whose  axes  are  allowed  to  vary  independently,  thus  modelling,  to  first 
order,  the  changes  of  bubble  shape.  The  Lagrangian  equations  of 
motion,  incorporating  the  effects  of  the  change  of  shape  and  of  energy 
loss  by  the  radiation  of  sound,  are  derived  and  solved  for  the  case  of 
a  spheroidal  bubble.  -The  results  of  these  calculations  for  various 
initial  conditions  are  compared  with  analgous  cases  for  a  spherical 
bubble. 

»  / 

It  is  found  that  the  spheroidal  bubble  model  predicts  a  reduction 
In  the  upward  translational  motion  of  the  bubble^of  a  factor  of  approx¬ 
imately  2.  A' comparison  of  the  predicted  upward  motion  of  a  spheroidal 
bubble  produced  by  227.27  kg.  of  TNT  detonated  46  metres  below  the 
surface  shows  good  agreement  with  that  which  has  actually  observed. 
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I.  INTRODUCTION 


The  formation  of  a  bubble  of  gaseous  detonation  products  always 
accompanies  an  underwater  explosion.  This  bubble  rises  toward  the 
surface  of  the  water,  responding  as  it  does  so  to  the  change  In  the 
external  pressure  distribution  with  oscillatory  motion,  during  the 
course  of  which  It  loses  some  of  Its  energy  through  the  emission  of 
sonic  pulses.  Although  the  bubble  Is  Initially  spherical,  the  effect 
of  Its  upward  motion  Is  to  distort  It  into  a  non-spherlcal  shape,  which 
becomes  most  pronounced  In  the  neighbourhood  of  the  minimum  radius. 

The  alteration  In  the  bubble's  shape  further  affects  both  the  pulsa¬ 
tions  and  the  upward  translational  motion  of  the  bubble.  By  means  of 
finite  element  techniques,  the  equations  of  motion  of  the  bubble  can  be 
solved,  taking  Into  account  the  effects  of  the  changing  shape  of  the 
bubble,  although  the  amount  of  computing  time  required  by  this  method 
limits  Its  utility.  Finite  element  methods  have  a  further  disadvantage 
In  that  physical  Insights  into  the  systems  considered  are  rather  more 
difficult  to  come  by  than  might  otherwise  have  been  the  case. 

Herring  (1942)  and  others  (eg.  Taylor  (1942)  and  Shiffman  and 
Friedman  (1944))  have  treated  the  problem  of  the  motion  of  the  bubble 
by  considering  It  to  be  a  perfect  sphere  throughout  its  entire  motion. 
This  treatment  yields  values  for  tlie  periods  of  radial  pulsations  of 
the  bubble  which  are  in  good  agreement  wltli  experimental  data,  but 
predicts  a  much  more  rapid  movement  toward  the  surface  than  Is  actually 
observed.  This  arises  because  the  largest  upward  velocities  of  the 
bubble  occur  at  tliose  times  when  the  bubble  is  near  its  minimum  radius; 
It  is  precisely  then  that  the  largest  departures  from  sphericity  occur. 
Penney  and  Price  (1942)  and  Ward  (1943)  included  the  effects  of  the 
non-sphericity  of  the  bubble  on  Its  motion.  However,  It  was  always 
explicitly  assumed  In  their  derivations  of  the  velocity  potential  of 
the  flow  about  the  bubble  that  the  departures  from  sphericity  are 
always  small.  Accordingly,  tlielr  equations  are  not  applicable  near  the 
times  at  which  the  bubble  is  at  Its  smallest  volumes. 

Hicks  (1972)  was  able  to  bring  the  value  for  the  upward  transla¬ 
tional  velocity  of  a  spherical  bubble  at  a  minimum  radius  into  agree¬ 
ment  with  experimental  data  by  adding  a  drag  term  to  tlie  equations  of 
motion.'  In  his  formulation,  the  drag  coefficient  Is  an  empirical 
correction  whose  value  Is  cliosen  to  make  the  predicted  rate  of  rise  a!; 
the  first  minimum  radius  consistent  with  observation.  However,  a 
different  drag  coefficient  must  be  selected  for  each  cliarge  mass  and 
depth,  requiring  a  comprehensive  data  base  from  which  the  appropriate 
value  can  be  chosen  for  each  case.  For  these  reasons,  a  model  for  the 
bubble  In  whlcli  large  deviations  from  sphericity  and  their  effects  or. 
Its  translational  motion  are  treated  would  be  of  considerable  practic-.l 
and  theoretical  interest. 


Ill  this  work,  a  Lagrani»ian  is  derived  for  a  bubble  whose  shape  is 
not  constrained  to  be  always  spherical  but  may  become  ellipsoidal  as  it 
moves  upwards.  The  equations  of  motion  for  a  general  ellipsoidal 
bubble,  incorporating  the  effects  of  loss  of  energy  by  radiation,  are 
presented.  The  algorithm  by  which  these  equations  of  motion  were 
solved  numerically  is  briefly  discussed.  Computational  results  for 
some  charge  masses  and  depths  are  presented  and  compared  with 
experimental  data. 

II.  REVIEW  OF  PREVIOUS  WORK 

Taylor  (1942)  derived  equations  describing  the  motion  of  a 
spherical  babble  of  gas  undergoing  both  radial  pulsations  and 
translational  motion  toward  the  water  surface.  These  are: 


2Tipa^  (^) 


[2.1] 


=  -  E(a) 


"  dT  -  7T  i  ^ 


f2.2] 


where  a  Is  the  radius  of  the  bubble  as  a  function  of  the  time  t,  U  its 
upward  velocity,  z  the  position  of  the  bubble  below  the  pressure  datum 
(i.e.  below  the  zero  pressure  level),  E(a)  the  Internal  energy  of  the 
gas  comprising  the  bubble,  g  the  gravitational  acceleration,  the 
total  energy  of  the  bubble  and  p  the  density  of  the  water.  In  Taylor’s 
formulation,  there  was  no  mechanism  Included  for  energy  loss,  and  hence 
Y^^  was  taken  to  be  a  constant.  For  TMT  explosions,  it  has  been  found 
(Herring  1942)  that  approximately  50%  of  the  total  explosion  energy  is 
retained  by  the  bubble;  in  that  case. 


=  (1.85  X  10^'^);i 


[2.3] 


where  Y^  is  measured  in  ergs,  and  M,  the  original  mass  of  the  explosive 
charge,  is  given  in  gra.  If  one  assumes  that  tlie  gaseous  explosion 
products  obey  the  ideal  gas  law,  then  the  internal  pressure,  P,  is 
given  by 


P  = 


[2.4] 


where  p  is  tiie  density  of  the  explosion  products  and  y  t'ne  ratio  of 
specific  heats.  Assuming  that  the  entire  mass,  M,  of  the  explosive  has 
been  converted  to  gas. 


and  hence, 


[2.6] 


[2.7] 

/  =  1.25 

for  TNT  where,  In  eq.  [2.6],  E(a)  is  measured  In  ergs  and  m  in  gm. 

Hicks  (1^72)  Lncoi-porated  a  drag  force,  into  the  equations  of 

motion,  where  ifj  given  hy 

~  T  CpjTtpa^U^  [2.8] 

and  the  value  of  the  drag  coefficient,  Cp,  was  chosen  to  he  =  2.25 
in  order  to  bring  tlie  distance  travelled  upward  by  the  bubble  at  its 
first  maximum  into  agree. iiet\c  with  that  actually  observed  for  500  lbs. 
of  TNT  detonated  150  ft.  below  the  surface.  By  differentiating  eq. 
[2.2],  he  obtained  the  rate  of  change  of  moraentura  wltli  respect  to  time. 
The  incorporation  of  Fj^  into  the  equation  of  motion  yielded: 

^a^U)  =  2a^g  -  i  C^a^U^  [2.9] 

L,  the  Lagrangian  of  the  bubble,  is  given  by 

L  =  2xpa^(-^)  +  j  oa^U^-  -  ^-^a^pgz  l-'l^O] 

-  E(a) 

Hence,  a  more  general  form  for  the  equations  of  motion  of  a 
spherical  bubble  is  given  by 

d  /  hi.  _ 

-0,,. 


CO 

E(a)  =  /  PdV 
<1 

.  kiiL€21r}2 
(y-1) 

•y 

where  dV  =  Aua'^da.  Taylor,  using  the  work  of  Tones,  set 
k  =  7.83  X  10" 


= 


[2.11] 

-U,  and  the  Qj^  are  the  generalised  dissipative 


where  a  = 
forces . 


The  (I  issiiJ.it  Ive  draii  force,  is  given  by  Fp  in  eq.  [2.8).  It 

has  prevloii.sly  been  shown  (Heaton  1984)  that  the  generalised  force 

for  dissipation  by  radiation  of  sound.  Is  given  by 


2  i2  _  J_ 


Qa  =  ^  -  t 


[2.12] 


where 


AQ  =  4iipn^  [  4.aaa*  +  (a“a*“/a)] 


[2.13] 


a  =  y,  and  Cg  is  the  speed  of  sound  in  the  water, 
dt 

The  terms  contained  in  AQ^  arc  analogous  to  the  radiation  reaction 
terms  In  electromagnetic  theory,  and  hence  can  be  Ignored  in  a  first 
approximation,  although  this  approach  will  undere.stlnate  the  energy 
loss  near  the  minimum  r.idius,  and  overestimate  It  elsewhere. 


Ward  (1943)  and  Penney  and  Price  (1942)  derived  equations  of 
motion  for  a  ne.arly  spherical  bubble  by  expanding  the  velocity  poten“ 
tlal,  <5,  of  Che  flow  about  the  bubble  in  terms  of  the  Legendre  poly¬ 
nomials,  tiiusly: 


’’KcosO)  ^  ^(cos9)  ^ 


[2.14] 


where  the  coefficients  A,  and  B2  are  functions  only  of  time.  They 
further  assumed  that  the  radius  vector,  R(t),  from  the  centre  of  the 
bubble  to  a  point  on  Its  surface  could  be  written  as 


R(t)  =  a  +  b2P2(cosG)  +  b3P.j(cose)  +  ... 


I  2.15  i 


where  a,  hj,  and  b'^  are  functions  of  time  only.  At  the  surface  of  the 


bubble , 


S  =  -  (I?)  -  " 


[2.16] 


=  ^  +  2  — 4  Pi(cosO)  +  3  — T-  Pt(cos0) 
R^  rJ  R 


Substituting  eq.  |2.15l  Into  eq.  [2.16]  and  equating  the  coefficients 
of  Che  Legendre  polynomials  on  both  sides  of  the  equation  yields; 


^2  da 

^  77* 


dbo  \ 

^^2^2  = 


[2.17] 


OUT  \  '-"a  I  g  . 
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the  surfneo  of  the  bubble,  tlie  pressure  must  be  iinffarm  and  e'|iial  to 
the  Internal  ^as  pressure.  ilsLn^j  Bernoulli's  equation,  this  condition 
can  be  written  as 

2 

gz  -  gRcos9  +  ~  i(  =  |(P8)^  [2-18] 

If  one  substitutes  for  ■!>  and  R  in  eq.  [2.18],  using  eqs.  [2.14]  - 
[2.15],  multiplies  the  resulting  equation  In  turn  by  each  of  the 
Legendre  polynomials,  and  then  integrates  over  cos0,  the  orthogonality 
relations  among  the  Legendre  polynomials  produce  4  differential  equa- 
t Ions : 


gz  +  a(^f)  ^  J  ^  +  0(uS  =  ^(PJY. 


dt^ 
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da 
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0 
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a 
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U“  - 
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a“  dt“ 
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da 

1 
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;T  TT  ~ 
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Tt 

0  7^ 

1 
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dt 


dt  +  OCU"*) 


[2.19] 


Ward  (1943)  has  estimated  the  value  of  82*  which  measvires  the  depar¬ 
tures  from  sphericity,  in  eq.  [2.15],  and  found  that  it  remains  small 
until  tlie  bubble  begins  to  contract.  Near  the  bubble's  mlnlraun  radius, 
b2  beco.nes  greater  than  a,  making  the  wliole  calculation  invalid. 

Ill.  EQUATIONS  OF  MOTION  FOR  AN  ELLIPSOIDAL  BUBBLE 


Now,  let  a,  b,  c  be  the  scini-axes  of  an  ellipsoid  along  the  x,  y, 
axes,  respectively,  for  a  co-ordinate  system  whose  origin  is  at  the 
centre  of  the  ellijjsoid.  Let  t!ic  ellipsoid  be  immersed  In  a  fluid  of 
infinite  extent,  and  let  one  of  the  axes,  say  a,  vary  witli  respect  to 
time.  At  any  Instant  in  time,  the  equation  of  the  ellipsoid  will  bo 
given  by 


jC 

a 


2 

2 


+ 


3.1 


The  velocity  iiotential,  T’,  for  the  flow  about  the  ellipsoid  is 
given  by  the  solution  to  Laplace's  e<iuation, 

V“T>  =0  [3.2] 

with  .ipproprlate  boundary  conditions.  Since  the  problem  obviously 
possesses  ellipsoidal  symmetry,  it  is  most  convenient  to  transform  to 
ellipsoidal  co-ordinates,  thusly: 
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[3.3] 


,2  .  iaimuit4(4+ii), 

(n^-b2)(,2.j2) 

2  .  [3.,] 

,2  ^  i£!±U(Ztui£c^ 

where  \,  n,  v,  are  Lhe  ellipsoidal  co-ordinates.  The  surfaces  defined 
by  X,  =  constant,  p  =  constant,  v  =  constant,  are  confocal  quadrics. 
Because  of  the  symmetry  which  exists  In  the  transformation  equations, 
eq.  [3.3],  one  is  allowed  to  specify  which  co-ordinate's  constancy  will 
yield  confocal  ellipsoids.  Throughout  this  paper,  then,  the  relation  \ 
=  constant  will  be  taken  to  describe  the  family  of  confocal  ellipsoids. 
In  ellipsoidal  co-ordinates ,  Laplace's  equation,  eq.  [3.2],  is  given 
by: 


,  2 
4-)  -h 


^  (v->0(k.^  -2^) 

where  k.^  Is  given  by 

1 

=  ((a^+X)(b2+\)(c“+X))'' 


[3.4] 


[3.5] 


and  the  expressions  for  and  k^  can  be  obtained  from  eq.  [3.5]  by 
symmetry. 

Now,  Let  a  he  a  solution  to  eq.  [3.4],  and  let  another  solution 
be  given  by 


[3.6] 


where  x  Is  ^  function  of  X  only.  By  substituting  eq.  [3.6]  into  eq. 
[3.4],  one  finds  that  a  must  have  the  form 

■'i  =  '!’>?(  fi.v)  [3.7 


where  is  a  function  of  X,  only  (Mllne-Thompson  1949).  Using  eq. 
[3.6]  to  aid  in  the  solution  of  eq.  [3.4],  one  finds  that 


X(X)  =  A  /  +  B 


[3.3] 


wliere  A  and  B  are  arbitrary  constants.  Hence,  if  a  is  a  solution  to 

eq.  [3.2], 

^  ^  a  j  [3.9; 

“x'^x 


is  also  a  solution.  The  solutions  to  eq.  [3.4]  having  the  Cona  of  eq. 
[3.9]  are  the  ellipsoidal  harmonics. 

In  the  case  of  the  spherical  bubble,  Taylor  and  Herring  assumed 
that  its  pulsations  would  be  described  by  simple  radial  oscillations. 
Since  a  sphere  is  a  degenerate  ellipsoid,  it  would  be  reasonable  to 
choose  as  a  solution  to  eq.  [3.4]  the  ellipsoidal  harmonic  which  pro¬ 
duces  analogous  oscillations.  Such  a  velocity  potential  is: 

«  -  A  I  ^  [3-10] 

since  ■  1  is  a  solution  to  eq.  [3.4].  The  upper  limit  of  the  inte¬ 
gral  has  been  clioscn  in  order  that  the  potential  become  0  at  an 
infinite  distance  from  its  source.  The  constant  of  Integration  A  is 
determined  by  the  boundary  conditions.  Now,  the  boundary  conditions 
for  a  pulsating  ellipsoid  are  not  as  obvious  as  those  for  a  sphere. 
Nonetheless,  it  seems  apparent  that  at  that  point  on  an  axis  which  is 
on  the  surface  of  the  ellilpsold,  the  normal  velocity  of  the  fluid  must 
be  equal  to  the  rate  of  change  with  respect  to  time  of  that  axis. 

Hence,  for  an  ellipsoid  in  whlchonly  one  axis  is  allowed  to  vary, 

say  a, 

©  -- 

x^a 

y-o 

z-o 

where  a  ■  ^  ••  ^$»n  is  the  normal  derivative  of  the  potential,  and 

n  is  tlie  unit  outward  normal  to  the  ellipsoid.  In  ellipsoidal 
coordinates,  eq.  [3.II]  can  be  written  as: 

2^^)  -a  [3.12] 

(pv)i 

x=a 

y-o 

z«o 

where  the  potential  $  is  given  by  eq.  [3.10]. 

How,  at  x»a,  y=o,  z"0 


[3.11) 


[3.13] 


and 

■  abc 


[3.14] 
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A  =  -y  abc 


[3.15J 


i  abc  I 


[3.16] 


Tlie  kinetic  energy  oF  the  fluid  around  the  ellipsoid  is  given  by 


^  P  #  ^ 


[3.17 


where  the  integral  is  carried  out  over  a  bounding  ellipsoid  described 
by  eq.  [3.l],  and  over  a  second  outer  boundary  which  is  obtained  by 
allowing  \  in  eq.  [3.3]  to  approach  Infinity,  essentially  describing  an 
iifinitely  extended  ellipsoid.  The  contribution  to  eq.  [3.17]  from  the 
large  ellipsoid  vanishes,  leaving  only  that  at  the  inner  surface  i.e. 
the  bubble  itself.  Tne  surface  integral  can  ho  transformed  into  one 
over  t!ie  x-y  plane,  thusly: 


^^■PdS  =  //  dx  dv 

^  3n  bn  (-.5) 


[3.18] 


where  z  is  the  unit  normal  along  the  z  axis.  The  unit  normal,  n,  to 


the  surface  of  tt>e  ellipsoid  is  given  by 

^  1  /XA.vA.zAs 

n  - - T  X  +  y  +  — s-  z) 

I  z  _  4  ‘  a  b  c 

+  ■“t) 

^a^  b"' 

in  the  ellipsoidal  co-ordin  ites , 

_  L 


[3.19' 


A.=  0  L  \=0 


whe  re 


1.2  _  i  /  X  ,  y“  .  z-  N 

^  ^  (a-+X)“  (b  +X)^  (c‘^+\)- 

So,  at  the  surface  of  the  ellipsoid,  \=o. 


1 

- 1 

2  2  2  ■ 


[3.20] 


'3.21 1 


,3.22] 


Using  eqs.  i3.16],  [3.18],  !3.19|,  and  [3.22]  and  taking  account  of  the 
contributions  from  t'ne  half  surfaces  above  and  below  the  x-y  plane,  eq. 
[3.17  ]  becomes 


Evaluating  eq.  1^3.23]  finally  yields  the  expression  for  the  kinetic 
energy,  T,  of  tlie  flow  about  an  ellipsoid,  described  r,y  eq-  [3.1  ],  wlien 
the  semi-axis  a  is  allowed  to  vary  with  resjiect  t'^  time: 


T  =  npa“b^c^  I ”4^  13.23] 

Since  the  equations  are  symmetrical  with  respect  to  all  three  semi-axes, 
it  is  possible  to  obtain  expressions  for  tlio  Icinetic  energy  of  the  flow 
about  an  ellipsoid  when  each  of  the  other  semi-axes  are  varying  b;'  means 
of  cyclic  permutations  of  a,  b,  and  c,  tiiusly: 


T  =  Tipab'^c*" 

r 

.4A 

[3.26] 

-T.  2.  2*?  ^ 

T  =  Tipa  b  c  ’  j 

0 

dX 

[3.27] 

whore  6  =  ^  and 


d>: 

dt'  * 


Hence,  by  adding  eqc.  [3.23 j  -  |3.27j,  and  renormalizing,  one  can 
obtain  the  kinetic  energy,  T,  for  the  flow  around  an  ellipsoid  when  all 
tiirec  axes  are  allowei  to  vary  in  time: 


= 


b"c‘ 


2^ 

a  n 


2,2 


dA 

k^ 


f3.28] 


The  normalisation  factor  3  has  been  chosen  so  that  eq.  [3.28]  is  in 
agreement  with  the  term  f(jr  the  kinetic  energy  due  to  radial  spherical 
pulsations  In  eq.  |2.l]  wiien  a=b=c  and  a=&=c. 


/\s  mentioned  above,  it  is  not  entirely  clear  what  sort  of  boundary 
conditions  are  applicable  at  tiae  surface  of  a  pulsating  ellipsoid.  in 
fact,  the  specification  of  boundary  conUitions  when  all  tiiree  semi-axes 
are  varying  is  equivalent  to  putting  constraints  on  the  inti^ract ions 
among  a,  fi,  and  c,  which,  in  turn,  is  equivalent  to  specifying  which 
stiapes  the  bubble  will  hi*  allowed  to  .issume.  Equation  [3.28]  corre¬ 
sponds  to  a  velocity  potential  wiiicli  has  been  so  constructed  that  the 
velocity  cross-terms  in  eq.  [3.17  cancel  out.  Since  a,  R,  and  c,  are 
all  mutually  perpendicular,  this  seems  pliysically  reasonable.  TViis  has 
the  effect  of  Insisting,  thal.  tlie  movement  of  a  point  on  an  axis  of  tlie 
ellipsoid  is  due  only  to  the  change  In  len'flh  of  t.iiat  axis,  with  the 
changes  in  tlie  lengths  of  the  other  two  axes  contributing  nothing. 


The  velocity  potential  associated  with  the  jxjrely  translational 
motion  of  the  ellipsoid  along  the  z  axis  is  well  known  (e.g.  Milne- 
Thompson  1949).  The  boundary  condition  is 


w  ) 

\=o 


Ucos0„ 


[3.29] 


where,  as  before,  -U  =  ~  ^2  angle  between  the  z  axis  and 

the  normal  to  the  surface  of  the  ellipsoid,  and  z  is  the  position  of  tlie 
bubble  centre  below  the  pressure  datum.  Since 


If  \,^ 

eq.  [3.29]  becomes 

\  -  ..  dz  \ 

\=o  "  '  \  =  o 


[3.30] 


’3-31  ] 


The  solution  to  eq.  [3.4]  which  satisfies  the  boundary  conditions 
in  eq .  [3.31  ]  is  the  first  order  ellipsoidal  harmonic  given  by 


$  =  Cz 


d\ 

Cc2+\)k, 


[3.321 


where  C  is  a  constant  of  integration.  Direct  substitution  of  eq.  [3.32] 
back  into  eq.  [3.31  ]  yields 


C  =  - 

'-«o 


'3-33] 


wliere 


a  =  a  be  /  — - 

°  o  (c^+X)k;^ 


[3.34] 


The  kinetic  energy  of  the  flow  around  a  translating  ellipsoid  is  then 
given  by: 


T  =  “TP/J'5-^dS 


f'(2~a")  //  cos0.,dS 


[3.35I 


=  abep 


As  before,  t!ie  Integration  is  carrietl  out  over  a  bounding  ellipsoid 
whose  semi-axes  are  a,  b,  and  c,  and  over  one  wViose  semi-axes  are 
allowed  to  extend  to  Infinity,  where  the  contribution  from  the  outer 
ellipsoid  vanishes.  Hence,  tiie  kinetic  energy  for  tlie  flow  generated  by 
a  translating  ellipsoid  whose  semf-axes  are  varying  in  tine  is  given 


o 


r3.36l 


.  2Tt  ,  f  ^o  vlj‘- 

+  -j  abcaCj--) 

o 

Equatiori  [3.3b  j  l\a;;  no  pravisioa  for  Lae  taLor.ict  i  oa  of  a,  fi,  or  c  with 
U.  Again,  this  is  not  unoipt'oto'i  siore  eq.  13.33]  should  yield  the 
kinetic  energv  terms  in  eq .  ‘2.1  !,  as  it  does  in  fact  do,  wlien  a=b=c  and 

;*  =  &=c  .  Phy:'.i  c.nl  ly ,  ns.  long  ns  the  bubble  is  sym::iet  r  ical  with  respect  to 
ttie  x-y  plane,  the  c.ont  rihut  ions  froi;  tlie  top  nnd  bottom  halves  to  a 
coupling  of  the  osc  i  1  l.ttory  motion  tjltli  ti;e  translat  iona  1  exactly  cancel 
each  other. 


The  energy,  V^,  n,:snciatec!  with  tlie  hydrostatic  pressure  around 
tlie  bubble  is: 

Vp  =  pgahcc  1 

The  internal  enei'gy  of  the  bubble,  r:(.n,b,c),  is  given  by 


h  ( a  ,  b  ,  c  )  - ; - -Tf - 

(y-1 


f  3 , 38  1 


Hence,  the  Lajrang  i  an,  L,  for  tlie  flow  around  an  ellipsoidal 
b'lliblo  whose  semi-axes  •■.V’'  i,  \i,  c,  at  n  dept.h  ^  lielow  tlie  pressure 
intum,  where  thrj  c  n<l'.  is  pirallel  to  tiie  ellipsoid  axis  c,  nnd  wiiich 
is  inoving  wit'n  <i  t  n  ns  1  a  t  ion  il  velocity  d  =  -  is  given  by 


L  =  ^  -(.•-t)"c'  h  a-lr’-’  ^  ^ 


^  TT  o  •  ^  *  A  ^ 

+  —  j  aabc  ~  — j  pabcgz 


f3. 39  j 


The  total  energy  of  r!ie  bubble,  Y(t),  at  any  time  t  is  therefore 
given  by: 

Y(t)  =  4  oCa^b^c^  .  a^b^c'^  -  a^h^c2) 

^  b  \ 


■h  — J  oabc(2-lr^)'  +  pabcgs 

(  Y~1  )  (“'"j) 


3.10  1 


Hence,  tlie  eqiiations  of  motion  for  Che  hubhle  .ire; 


'  i , \:t ■>  I  i  t ) 

n 


'>  >i^L  2i2(^L 

+  {2ab"l^  +  a  b  ^j^Jac  b  +  "*  ’  ^^)6c 

+  f-  0(<;2b2c2  b  b  a2l.2c-2) 

2 

+  p(a^b^c  +  pabt(a^)z 

2_  4  IT 

+  -j  pabc  -  z - j  pgabz 


.V_  /a.b.Cx*  2 

*  %a-a^ 

~  TT^ry 


wha  re 


•r  _  <l^a  V  _  d^b 
a - -T,  b  -  — T, 


[3.44 


[3.45] 


V  =  d^z 


[3.46] 


T  -  ^  d". 

■  i 


f(^n)  - 


-  rr-a-T* 


^:o  p  ^  .]£o  ? 

31  da  ^  5b  ^  ^  5c  ‘  » 


df(g^)  5t-(g^).  .  llL^o). 

dt  5a  5b  ^  5c 

.  =  kM Y _ 

In  order  to  complete  the  derivation  of  the  erjiiations  of  racjtion  of 
the  bubble,  It  Is  necessary  to  deternine  ()  ,  Q. ,  Q,  and  Q  ,  the  general- 

*3  o  z 

Ised  dissipative  forces  associated  with  the  radiation  of  sound  by  the 
bubble.  Now,  at  distances  large  compared  with  the  scale  of  the  bubble, 
the  form  of  the  velocity  potential,  'I’,  in  the  fluid  will  be  identical 
with  tliat  of  a  spherical  bubble.  Hence, 


'I-  =  -  ^  +  X  •  ?(l/r) 


[3.47] 


•  .■  -j.'- h'-  -V 


y.  .  y.y  ■  ,  .  n.  1^1, 


T*'«T*T*T»Tr' 


?  ■  'J  ■  11' 


wh(.!re  C  and  it  are  constants  (Vliioh  depend  only  on  the  time  t.  r  Is  the 
distance  to  tlie  field  point  frrxn  an  origin  located  somewhere  within  the 
bubble. 


Following  the  development  given  by  Landau  and  Llfshltz  (1966),  In 
the  wave  zone, 

[3.48] 

wliere  0  Is  the  angle  between  the  direction  of  the  translational  motion 
and  r,  and  the  retarded  time  t’  Is  given  by 


where  Cg  Is,  as  before,  the  velocity  of  sound  in  water.  The  velocity, 
iT,  of  the  water  In  the  wave  zone  must  therefore  be  given  by: 


[3.49] 


$ 


r^>  +  ^.(A(f)  ?) 


V  =  -  7$ 


c-r 


[3.5U] 


+  .  •  • 

where  terras  of  higher  negatWe  order  in  r  hiive  been  neglected.  The 

1  p 

total  energy  omitted  as  sonic  radiation  per  unit  time,  — -  is  then: 

dt 


^  -  pc,  //  (^•^)dS 

dt  b 

=  _  _  4Jia 

Cg  ht  3^3  ^^2 

9 


[3.51] 


where  the  integral  has  been  taken  over  a  sphere  of  radius  r  (Landau  and 
Llfshltz  1966). 

_o  To  a  good  approximation,  the  tenn  in  eq.  [3.51  ]  proportional  to 
C”^  can  be  neglected  for  low  translat ional  velijcities,  since  It  will  be 
S  —1 

2  orders  of  magnitude  smaller  tlian  that  proportional  to  Cg 


"low,  the  volume  4nC  of  fluid  which  flows  through  the  surface  over  which 
the  integral  in  eq.  f3.5l]  is  taken  must  be  equal  to  the  rate  of  change 
with  respect  to  time  of  the  volume,  of  the  bubble.  Thus, 


C 


_  1 


(abc  +  al*c 


+  a  he) 


[3.52] 


367 


-e.j  J’ 


aiu!  so 


[a be  +  abe  +  abc 
+  :i(  +  be) 

+  t(  ac  +  ac)  [3.53] 

2 

+  c( ab  +  afi) ] 


Now,  evidently, 

^  =  aQg  +  E.Q,^  +  ^  [3.54] 

Since  all  of  the 
was  contained  in 

same  approxinat ion ,  =0  in  eq.  [3.54].  Since,  in  the  absence  of  any 

translational  motion,  tiiere  exists  notiiing  to  distinguish  one  axis  of 
the  ellipsoid  from  another,  it  follows  that  one  should  be  able  to  obtain 
the  other  two  from  one  by  cyclic  permutation  of  the  axes  a,  b,  and  c. 
The  only  grouping  of  tiie  terns  in  eq.  ^3.53]  which  is  invariant  under 
cyclic  permutation  of  the  axes  is  given  by 


dependence  on  the  translational  velocity  in  eq.  [3.51 1 
the  term  proportional  to  Cc  it  follows  that,  to  the 


where 


F(a^) 


=  a.aja^  +  a^aja^ 
H-  a.  a..a  ,,  i^j^k 

1  N  j  ’ 


[3.56^ 


and  a,  =  a,  32  =  b,  a-j  =  c. 

Hence,  by  comparison  of  eq.  ^3.55]  with  eq .  [  3 . 54  j , 

[3.57] 

w'nere,  as  before,  i,  j,  k,  are  successively  equal  to  1 ,  2,  3.  Now,  as 
was  the  case  for  the  spherical  bubble,  the  terms  in  eq.  [3.57]  which 
depend  upon  tlie  products  ot  the  [>ulsational  accelerations  with  tliera- 
selves  or  with  the  pvil  sati^anal  velocities  are  analogous  to  the  radiation 
reaction  terms  in  alrc  tr  o  nay.net  Ic  theoiv.  This  suggests  that  such  terms 


i,  ^ 

i  •■> 


F(a.) 


.)-"i 


F(a  :  ) 


may  be  ij'iiorab le ,  at  least  ia  i  first  approximation-  At  this  stage,  in 
the  absence  of  any  perturbing  force  in  the  x-y  plane,  it  is  possible  to 
allow  6  =  a.  This  allows  one  to  drop  eq.  [3.42]  as  an  equation  of 
motion,  and  replace  R  and  h  by  a  and  a  in  eqs.  [3.43]  -  [3.45]  and  eq. 
[3.57].  This,  of  coarse,  specialises  the  equations  of  motion  to  those 
of  a  spheroidal  bubble.  Henceforth,  throughout  this  paper,  the  case  of 
the  spheroidal  bubble  will  he  treated  exclusively. 

To  sum  up:  e<js.  [3.42]  -  [3.45]  are  equations  of  motion  de¬ 
scribing  ii  pulsating  ellipsoidal  bubble  undergoing  translational  motion. 
When  Qa  =  Qb  “  Qz  “  equations  neglect  any  sort  of  energy  loss. 

It  has  been  shown  that  the  energy  loss  from  the  translational  motion  of 
the  bubble  can  be  expected  to  be  negligible  with  respect  to  that  from 
the  pulsational  motion  and  so  was  set  to  0  in  eq.  [3.45  j.  It  was 
further  shown  that 


p  3.1 


where 


F(a[)  =  a.Aja^  + 

+  AF(a»), 

When  one  wishes  to  ignore  the  effects  of  radiation  reaction, 
AF(a*)  =  0 


[3.58] 


L3-59] 


[3. 60] 


AF(a|)  ^  a[aja^.  istjT^k 
when  one  wishes  to  include  them. 


[3.61  1 


IV.  NUMERICAL  METHODS  OF  SOLUTION 

Before  one  ar.Lempt.s  numerical  solutions  of  the  equations  of 
motion,  eqs.  [3.43]  -  f3.45],  it  is  useful  to  make  them 
non-dimensional.  Thus,  tlie  subtst itut ion  of 


a  —  a  L , 
* 

c  =  c  L, 


[4.1] 


z  ^  L , 


t  =  t  T 


into  Lhe  equations  of  motion  used,  wliere 


L  =  (— ) 


|4.2] 


T=Vt- 


yields  a  dimensionless  form  of  the  equations  of  motion.  As  before, 
is  given  by  eq.  [2.3]  and  g  and  p  are,  respectively,  tlie  gravitational 
acceleration  and  the  density  of  water.  These  particular  scaling  factors 
in  eq.  [4.2]  were  originally  used  by  Taylor  (1942). 

Since  all  of  the  equations  of  motion  have  the  unfortunate  property 
of  singularity  at  the  origin,  it  is  necessary  to  begin  the  integration 
with  a  series  solution.  Taylor  (1942)  suggested  that  the  initial  solu¬ 
tions  to  the  dimensionless  forms  of  eqs.  [2.1  ]  -  [2.2]  be 


[^-3] 

*  5  t*2 

2  =  "'o  - 

where  z*  is  the^initial  dimensionless  deptli  below  the  pressure  datum, 
for  valc.es  of  t  near  zero.  Since  the  bubble  can^be  expected  to  be 
spherical  initially,  the  values  Cor  a  ,^z  ,  and  z  from  eq.  [4.3]  were 
used  to  beg^n  th^  integration  at  time  t^,  with  the  additional  require¬ 
ment  tiiat^a  =  c^.  One  also  needs  initial  values  for  tlie  rates  of 
change,  a  and  c',  of  t!ie  semi-axes.  These  were  estimated  by  assuming 
that  the  bubble  would  be  initially  spherlj^al  and  substituting  into  the 
d iraen| lonlesj  form  of  eq.  [2.1  ]  to  find  a  .  Hence,  the  initial  values 
for  a  and  c  are  given  by 


=  (1- 


fz'  .  2  * 


[4.4] 


•  •k  •  r 

a  =  C 


wliere 


>;*(..*) - viV.'25ri> - ,  ;4.5i 

(Y-l)(4n)Y  1 

Yg  is  the  initial  total  energy,  as  given  by  eq.  [2.3],  and  all  other 
variables  are  as  previously  defined. 

Another  numerical  difficiltv’  concerns  the  evaluation  of  the  terms 

^^o  ^“o  '^«o 

re"'  evidently. 


’ ««,-  V  •C' 


=  -  / 


fta 

SI  ® 

o  _  „  f 

- ^  J 

da, 


(1\ 

Trr  - 


hk  7  9  2  7  2 

^  (a-'+X)  (b^+X)  (c^+X) 

dX 

T - -  ^-r 


°  (a^+X)*  (b^+X)^  (c^+X) 
ail 


T77< 


-5^  =  bcl^  +  abc  -5^. 

9“o 


where 


[4.6] 

[4.7] 


dX 

1 - 5 - 2^ - ^ - T77 

(a-^+X)  (b^+X)  (c-^+X) 


[4.8] 


Iq  and  were  evaluated  with  the  TMSL  double  precision  subroutines 
MMLINF  and  MMLIND,  which  compute  incomplete  elliptic  integrals  of  the 
first  and  second  kind,  respectively.  The  partial  derivatives  of  1^,  eq . 
[4.6],  are,  in  fact,  incomplete  elliptic  integrals  of  the  second  kind, 
and  can  be  evaluated  with  MMLIND.  The  evaluation  of  the  second  terms  in 
eq.  [4.7]  presented  considerable  difficulty.  In  point  of  fact,  no 
commercial  routine  capable  of  evaluating  eqs.  [4.7]  seems  to  exist,  and 
the  difficulties  Involved  in  the  composition  of  one  ab  nihilo  are 
formidable.  Accordingly,  as  a  stopgap,  the  terms  in  eqs.  [4.7]  which 
involve  the  derivative  of  an  incomplete  elliptic  integral  of  the  second 
kind  were  evaluated  by  holding  one  of  a  or  c  constant,  and  varying  the 
other  at  each  step  in  the  integration,  tliusly: 

■^3“  °  TKa  -  (I]_(a-Aa,c)) ,  ['■'•9] 

ai,  , 

“  ZSc  c+Ac)  -  lj^(a,  c-Ac)) 


where  MMLIiJD  was  used  to  evaluate  Ij^.  Since  Aa  and  Ac  can  be  made  as 
small  as  desired,  theoretically  eq.  [4.9 J  can  be  made  to  approximate  the 
true  value  of  the  derivative  as  closely  as  desired;  however,  the  prac¬ 
tical  constraints  of  computational  time,  machine  accuracy,  and  the 
accuracy  of  the  IMSL  subroutine;:  do  place  limits  on  the  size  of  Aa  and 
Ac . 


The  actual  integrations  were  carried  out  using  a  4  point  Runge- 
Kutta  algorithm  incorporating  automatic  error  controls.  Some  numerical 
difficulties  with  this  method  were  encountered  when  the  radiation  re¬ 
action  was  incorporated  using  eq.  [3.61  ].  Near  the  minima  of  a  and  c, 
and  c  become  small,  and  change  sign  as  well.  Because  of  their  iepend- 
ence  upon  a~^  and  c 

adversely  affecting  the  c.onvergence  of  the  integration 
was  circumvented  by  the  use  of  a  series  approximation  in  which  a  ^  and 
c~^  were  replaced  by  averaged  pulsationnl  veloc.lties.  By  using  eq. 
[3.40],  it  is  possible  to  write 


,  the  values  of  and  can  oscillate  rapidly. 


This  difficultv 
-1 


•IT? 


jnJL-La-ia 

Y.,  -'ll 


[4. 10 


o  ’'u’'  ’'/"O  V  •*"  /  **  *  * 

-2a  b  c  (j^—)z  -  4a  h  c  z 


where 


*4  1  A 

a  =  ■gr  (a+b+c) 


-nrz  =  -1 
a  —I 


(a  “b  ‘■c  ^  +  a  -b  ‘• 


[4.11] 


+  a  b  c  ) 

and  Y(t)  Is  the  energy’  of  the  bubble  at  any  time  t,  given  by  eq. 

>■3.401. 

If  one  defines: 


f4.12i 


3  =  - 


2a*b*c* 


(2-^)2 


.2  /  *.*  *  * 

*  4  a  h  c  ^ 


then, 


■•*2)^  _  11 


^  ^  ^  ^  f.  ^  t...) 

where  the  positive  value  is  taken  while  a  particular  axis  is  expanding 
and  the  negative  while  it  is  contracting.  By  expanding  the  dimension¬ 
less  form  of  eq.  [3.58]  when  AF(a. )  is  given  by  eq.  [3.61  ],  and  sub- 

_  i  *• 

stituting  f a*^  1  for  (a*)  and  (c*)  values  for  the  dissipative 
function  incorporating  averaged  radiation  reaction  terras  were  obtained. 


Estimates  for  a  and  c  were  arrived  at  by  the  substitution  of  tlie 
current  values  of  a,  c,  a,  c,  z,  and  z  into  eqs.  [3.43]  and  [3.44]  with 
Q„  and  Q„  set  to  zero.  Those  values  for  a  and  c  were  substituted  back 
into  eq.  13.58|,  to  obtain  new  values  for  Q^and  which  were  in  turn 
used  in  eqs.  >^3.43]  and  [3.44]  to  obtain  new  estimates  for  a  and  c. 

V.  NUMERICAL  RESULTS  AND  ANALYSIS 


Figures  2-11  sliow  the  results  of  computations  for  a  bubble 
produced  by  the  detonation  of  2.1136  kg.  of  T'NT  6.1  metres  below  the 
surface,  using  eqs.  >3.43]  -  [3.46]  for  a  spheroidal  bubble.  Those 
curves  assoclate>d  with  a  spheroidal  bubble  and  labelled  'AY  =  O'  were 
calculated  tinder  the  assumption  of  no  energy  loss;  that  Is, 

=  0  in  eqs.  [3.42]  -  [3.46].  Tlie  curves  labelled  'AY  *  0,  AI'  -•  O',  were 


c  .lie  11  la  ted  laooriiorat  Lng  radlatlv/e  energy  loss,  but  not  radiation 
reaction  terns;  that  Is,  eqs.  [3.58]  -  |3.60]  wore  used  to  define  and 
Q^.  In  Figs.  2-5,  a  and  c  are  the  semi-axes  of  the  spheroidal  bubble, 
where  is  the  serai-axis  in  the  plane  nornal  to  the  bubble's  upward 

motion,  and  c  the  se.ni-axls  in  the  plane  parallel  to  the  bubble's  upward 
motion. 

Taylor  (1942)  considered  the  same  case,  using  eqs.  [2.1  ]  -  [2.2 j 
for  a  spherical  bubble  in  the  absence  of  any  energy  loss  from  any 
source.  In  Figs.  2-11,  the  curves  labelled  with  ' AY  =  O'  were  obtained 
by  solving  eqs.  [2.10]  -  [2.11  ],  with  “  Q2  “  which  are  equivalent 
to  Taylor's  equations  for  a  spherical  bubble.  The  curves  labelled  'AY  ^ 
0,  AQ^  =  0’  result  from  the  Incorporation  of  radiative  energy  loss, 
neglecting  radiation  reaction,  into  eqs.  [2.I0]  -  [2.11  ];  tiiat  is; 
was  given  by  eq.  [2.12],  with  AQ^  =  0.  The  curves  labelled  'AY  *  0, 

AQ.j  *  O'  i riciir p.oratG  energy  loss  Including  radiation  reaction;  that  Is; 
Q.^  and  AQg  wer.i  given  by  eqs.  [2.12]  -  [2.13].  In  Figs.  2-5,  'radius' 
refers  to  the  sidierical  bubble  radius  as  calculated  in  eqs.  [2.101  - 
'2.11]. 

The  curves  labelled  'spherical'  in  Figs.  6-11  are  the  upward 
velocities  and  heights  above  the  original  detonation  point,  obtained  by 
solving  eqs.  [2.10]  -  [2.11 J  for  a  spherical  bubble,  under  different 
assumptions  about  the  nature  of  the  energy  loss.  The  curves  labelled 
'spheroidal'  are  the  same  quantities  obtained  from  the  solution  of  eqs. 

I  3.43]  -  [3.45]  for  a  spheroidal  bubble. 

Taylor  (1943)  presented  photographs  showing  the  behaviour  of 
bubbles  generated  by  electrical  discharges  in  oil,  which  are  here 
reproduced  in  Fi,gs.  12-13.  Tliese  show  that  a  bubble  in  the  early  stages 
of  its  motion  is  very  nearly  spherical,  but  that  near  its  minimum 
volume,  it  becomes  approximately  disc  shaped,  with  its  longest  dimension 
lying  in  the  plane  normal  to  the  direction  of  its  upv;ard  motion.  Mear 
the  second  maxinim,  the  bubble  is  highly  non-spherical,  and,  in  fact, 
seems  to  be  atterajitin;;  to  fission,  exhibiting  an  extremely  large  bulge 
<?n  its  upper  surf, ice  and  a  flat  lower  surface.  After  the  bubble  has 
passed  through  its  second  maximum,  it  becomes  .mushroom-shaped  and 
actually  does  bifurcate  -at  its  second  minimum.  The  two  halves  rejoin 
later  to  form  a  distorted  disc. 

As  one  can  see  from  Fig.  2,  a  spheroidal  bubble  reproduces  the 
salient  features  of  the  observed  behaviour,  at  least  qualitatively.  a 
and  c  were  very  nearly  equal  to  each  otlier,  as  well  as  to  the  spherical 
radius,  at  the  first  maximum.  Near  the  first  minimum,  tlie  bubble  became 
more  obviously  spheroidal,  with  the  ratio  a/c  .assuming  a  value  of  2.34. 
After  the  first  minimum,  the  qualitative  agreement  between  the  sphe¬ 
roidal  bubble  model  and  Taylor's  photographs  was  less  pronounced.  This 
was  hardly  surprising,  giviui  that  1  he  equations  of  motion  constrain  th.e 


possible  sli.ipes  of  Llie  bubble  to  spheroids.  however,  the  model  dies  .it. 
le-ist  predict  that  the  bubble  would  not  retur.i  to  a  spherical  stMp.j.  In 
addition,  the  splieroidal  model  also  ret.ains  one  of  the  successful 
features  of  the  spheric. il  model,  that  of  the  prediction  of  the  period  of 
the  bubble's  oscillation.  The  period  of  the  fir.st  oscillation  in  Fig.  2 
(wliich  was  defined  for  both  the  spherical  -nnd  spheroidal  bubble  to  be 
the  time  after  a  maximum  at  which  the  first  derivative  of  the  volume 
changed  sign)  was  .26  seconds  for  both  models.  It  should  also  be  noted 
that  the  time  at  whicii  the  volume  was  a  minimum  did  not  coincide  with 
the  minima  of  either  a  or  c.  In  Fig.  2,  a  reached  Its  first  minimum  a 
full  .01  seconds  before  c  did. 

Figures  3-4  show  the  effects  of  Incorporat  inp,  rcdlatlvo  energy 
loss  into  the  spheroidal  bubble  model,  and  compare  the  result';  with  the 
analogous  case  for  a  spherlc.al  bubble.  When  radiative  energy  loss  w.is 
included  In  the  c.nlculat  1  ons ,  the  second  minima  of  a  and  c  occu-red 
slightly  earlier,  and  was  more  nearly  coincident.  Figure  5  compares  trie 
semi-axes  for  a  spheroidal  bubble  under  various  assumptions  abouf  Llie 
form  of  the  energy  loss.  It  Is  Interesting  rh.at  the  most  noticeable 
difference  .among  the  calculations  was  In  t'no  val'ie  of  c. 

The  outstanding  f.allure  'jf  the  si)herical  model  for  the  bubble  is 
its  prediction  of  .a  too  rapid  r.ate  of  rise  when  the  bubble’s  volume  is  a 
minimum.  As  Figs.  6-10  show,  the  m.i-timum  upw.ard  velocity  predii:ted  by 
the  spheroidal  model  was  less  than  Liiat  predicted  by  the  spherical  model 
by  a  factor  of  2,  and  consequently  the  distance  travelled  from  the  site 
of  Che  explosion  was  decreased  by  about  the  s.ime  amount.  This  diminun- 
t  ion  of  the  upward  translational  velocity  of  the  bubble  near  Its  mi  ninuTi 
volume  accounted  for  the  dif f erenci.'S  between  the  periods  of  tlie 

r !  cal  and  siiheroldal  models  -after  the  fir.st  minimum  in  Figs.  2-4. 
Recausi:  the  sj)liero  idal  buhljle  was  dci'per  than  the  spherical  one,  Llie 
hydrostatic  pressure  w.is  greater,  making,  the  period  shorter.  This  is 
also  Che  reason  Chat  Che  curves  I'l  ilgs.  2-11  associated  witli  the 
spherical  bubble  model  terminated  .1  second.s  before  those  of  the  sphe¬ 
roidal  model.  Because  of  t’ae  spheric.al  bubble's  greater  upward  veloc¬ 
ities,  it  readied  the  surface-  bef-irc  the  spheroidal  one. 

Figure  11  .sliows  Che  energy  possessed  by  tlie  bubble  as  a  fuuctioa 
of  time  for  botli  the  spherical  and  splier-jidal  models.  Since  tin-  iriir  ial 
solutions  to  the  equations  of  motion  for  the  spheroidal  bubble  were 
calculated  by  assigning  it  to  li.ive  been  initially  spherical,  it  is  not 
surprising  that,  !n  the  absence  of  dissipation,  tlci  energies  predicted 
by  the  two  models  were  found  to  be  constant  .uiii  eipial.  The  energv 
losses  predicted  by  tiie  spiu.'roid.al  model  without  r.idiation  re-acti'a  w.-re 
f'lund  to  be  in  close  .ig,reement  with  tlios.i  of  llie  spherical  model  without 
ridlatlon  reaction.  After  the  first  minimum  volume  was  passed,  Che 
spheroidal  model  predicted  a  slightly  gre.ater  energ.y  loss  than  the 
spherical,  until  the  spheroidal  bubble  p.i'ised  through  a  second  minimum 
volume . 


Tlie  most  dramatic  diifiironce  between  the  spheroidal  and  spherical 
models  is  In  the  role  of  radiation  reaction.  As  can  be  seen  in  Fi,>.  11, 
tlie  inclusion  of  radiation  reaction  in  the  splierical  bubble  model 
decreased  the  amount  of  energy  radiated.  For  the  spheroidal  bubble, 
this  was  the  case  only  until  the  first  minimum  volume  was  passed.  After 
that  time,  the  effect  of  including  radiation  reaction  In  the  calculation 
was  to  increase  the  radiative  energy  loss,  compared  both  to  that  of  the 
spheroidal  model  without  radiation  reaction  and  to  that  of  either 
spherical  model.  Moreover,  since  even  without  the  inclusion  of  radia¬ 
tion  reaction,  the  spheroidal  model  yielded  a  greater  energy  loss  than 
the  analogous  spherical  case,  it  seems  that  this  was  not  an  artefact  of 
the  approximations  used  in  the  computation  of  the  radiation  reaction 
terms.  lihether  the  magnitude  of  the  Increase  in  the  energy  loss  which 
occurred  when  radiation  reaction  terras  v^;ere  added  to  the  equations  would 
be  as  large  as  that  indicated  by  Fig.  11  is  rather  more  uncertain.  The 
close  correspondence  between  the  predicted  energy  losses  for  the 
spherical  and  spheroidal  models  prior  to  the  first  minimum  suggests  that 
the  calculation  was  valid,  at  least  in  the  regime  in  which  the  bubble 
was  nearly  splierical.  However,  as  mentioned  above,  after  the  first 
minimum,  the  actual  shape  of  the  bubble  is  not  really  spheroidal. 
Consequently,  the  actual  radiation  luss  may  be  quite  different  to  that 
calculated  for  a  spheroidal  bubble.  However,  since  the  bubble  is  even 
less  spherical  than  it  is  spheroidal,  on  balance,  it  seems  probable  that 
the  predictions  of  the  spheroidal  model  were  more  accurate  than  those  of 
the  spherical. 

Hicks  (1972)  solved  the  equations  of  motion  for  a  spherical 
bubble,  with  the  addition  of  a  hydr<}dynamlc  drag  term,  for  227.27  kg.  of 
TMT  at  a  depth  of  45.73  metres  below  the  surface.  (That  is,  Hicks  tot k 
eqs.  [2.11]  as  his  etjuatlons  of  motion,  witli  as  defined  by  eq.  |2.8j 
and  Qjj  =  0)-  Tt  has  been  observed  tl\at  a  bubble  from  an  explosion  with 
these  characteristics  rises  approximately  3.35  metres  from  the  location 
of  the  explosion  in  the  time  taken  to  reach  its  first  minimum.  Hicks 
found  that  where  drag  is  the  only  source  of  dissipation,  a  drag  coef¬ 
ficient  of  Gq  =  2.25  had  to  be  introduce!  into  the  equations  of  motion 
for  a  spherical  bubble  in  order  to  reproduce  tliis  behaviour.  In  this 
work,  for  tlie  initial  solutions  chosen,  it  was  found  that  a  drag  coef¬ 
ficient  of  Cp  =  i.85  brought  the  predicted  rise  of  a  spherical  bubble 
into  better  agreement  with  observation.  When  radiative  dissipation  was 

also  Included,  a  drag  coefficient  of  Cj^  =  1.6  yielded  better  agreeement 

with  observation. 

In  Figs.  14-23,  the  curves  labelled  'spherical'  were  olitaLneJ  by 
taking  eqs.  [2.IO]  -  [Z.ll]  as  the  cijuations  of  motion  for  a  spherical 
bubble  produced  by  227.27  kg.  of  TNT  detonated  45.73  metres  below  the 
surface.  All  of  the  symbols  in  the  legend  for  these  figures  have  the 
same  meaning  as  in  Figs.  2-11.  It  should  also  be  understood  that,  for 
the  spherical  bubble,  the  effects  of  drag  were  ignored  unless  a  valiie  of 
the  drag  coefficient  is  given  in  the  legend.  'When  lirag  was  con¬ 
sidered,  Ff^  was  given  by  eq.  '2.8].  Tt  should  bo  pm|)iias i zed  that  a  drag 
terra  was  Incorporated  only  Into  the  equations  of  motion  for  a  spherical 

bubble  and  never  Into  eqs.  r3.43l  -  [3.46],  the  equations  of  motion  for 

n  spheroidal  bubble. 


The  behaviour  of  the  spheroidal  bubble  in  Figs.  14-17  was  very 
similar  to  that  in  Figs.  2-5.  One  difference  was  tlie  close  agreement 
among  all  four  curves  for  the  period  of  the  first  bubble  oscillation. 

As  well,  for  tlie  second  oscillation,  the  period  predicted  by  the  sphe¬ 
roidal  model  was  closer  to  tliat  of  the  spherical  model  wlLli  drag  than 
that  which  was  predicted  by  the  spherical  model  without  drag.  This 
applied  in  the  cases  for  which  radiative  energy  loss  was  considered  as 
well  as  those  cases  for  which  it  was  not.  In  Figs.  18-21,  it  can  be 
seen  that  the  peak  velocities  predicted  by  the  spheroidal  models,  both 
with  and  without  radiative  energy  loss,  were  in  agreement  with  those 
predicted  by  the  spherical  model  with  drag,  most  notably  at  the  first 
minimum. 

It  is  Fig.  22,  liowever,  which  demonstrates  the  accuracy  of  the 
spheroidal  model.  As  can  be  seen,  the  height  above  the  site  of  the 
original  explosion  predicted  by  the  spheroidal  model  was  in  good  agree¬ 
ment  with  that  predicted  by  the  spherical  model  with  drag,  while  that 
predicted  by  the  spherical  model  in  the  absence  of  drag  disagreed  with 
that  predicted  by  the  spherical  model  with  drag.  Because  the  drag 
coefficient  used  with  the  spherical  model  was  chosen  specifically  to 
force  the  calculated  height  to  agree  with  observation,  the  agreement  of 
the  spheroidal  model  with  the  spherical  model  in  this  case  consitutes  a 
verification  of  the  spheroidal  model.  The  agreement  was  not  as  good  at 
the  second  minimum,  especially  for  the  spheroidal  bubble  without  energy 
loss.  However,  once  radiative  energy  loss  was  added  to  the  spheroidal 
model,  the  curves  exhibited  close  agreement  with  those  for  the  spherical 
model  with  drag. 

Figure  23  compares  the  predicted  energy  losses  of  the  spherical 
and  spheroidal  models.  Once  again,  it  can  be  seen  that  the  inclusion  of 
radiation  reaction  terms  increased  the  predicted  energy  loss  for  a 
spheroidal  bubble,  in  contrast  to  the  spherical  model,  for  which  the 
addition  of  reaction  terras  decreased  the  energy  loss.  It  is  inter¬ 
esting,  however,  the  energy  remaining  to  the  spheroidal  bubble  was 
greater  tlian  that  remaining  to  the  spherical  bubble  wlien  losses  from 
both  drag  and  radiation  were  Included. 

VI.  CONCLUDING  REMARKS 

In  this  work,  a  Lagranglan  for  an  oscillating  ellipsoidal  bubble 
which  is  also  undergoing  translational  motion  has  been  derived,  and 
equations  of  motion  obtained  from  it.  An  expression  for  the  generalised 
dissipative  forces  caused  by  the  radiation  of  sound  by  the  bubble  was 
also  found,  and  Incorporated  into  the  equations  of  motion.  This 
equation  was  specialised  to  the  case  of  a  spheroidal  bubble,  and  the 
equations  of  motion  solved  for  some  different  charge  m.asses  and  depths, 
both  with  and  wltliout  the  effects  of  radiation  of  sound  having  been 
Included . 


By  comparison  wltli  tiie  results  obtained  from  Taylor's  splierical 
bubble  model,  It  was  shown  that  the  spheroidal  model  retained  tlie 
successful  features  of  the  spherical  model,  notably  the  prediction  of 
the  bubble's  oscillatory  period,  and  reproduced,  at  least  qualltat Wely , 
notable  features  of  a  real  bubble's  behaviour,  such  Its  flattening  near 
Its  first  minimum,  and  a  slower  rise  time  than  that  predicted  by  tlie 
spherical  model.  After  the  first  minimum,  the  bubble's  shape  was  not  as 
well  modelled  by  a  spheroid;  however,  the  results  obtained  at  those 
times  were  still  superior  to  those  from  the  splierical  model,  particu¬ 
larly  with  respect  to  the  rise  time.  For  tlie  case  cited  by  Hicks 
(1972),  the  spheroidal  bubble  model  produced  results  which  were  in  very 
close  agreement  wltli  experimental  data.  The  spherical  model  was  capable 
of  similar  agreement  only  for  a  limited  t  line  and  only  with  the  addition 
of  a  drag  term.  The  spheroidal  model's  advantage  arises  from  its 
reproduction  of  the  height  above  the  explosion  site  naturally,  without 
tlie  addition  of  a  drag  term  wlilch  must  be  determined  from  experimental 
data  for  each  case.  Incidentally,  because  the  spheroidal  model  ignored 
liyd rodynamlc  dr;ig  completely,  and  still  reproduced  the  observed  be- 
liavlour  of  the  bulible,  it  seems  likely  that  drag  Is  relatively  unimpor¬ 
tant  In  determining  the  bubble's  motion;  It  appears,  rather,  that  the 
sha;x}  of  the  bubble  is  the  single  most  Important  factor. 

It  appears  that  the  spheroidal  model  predicted  a  gre;iter  loss  of 
energy  in  Llie  bubble  tlirough  tiie  radiation  of  sound  than  tlie  spherlciil 
model.  This  was  found  '"o  lie  most  significant  when  radiation  reaction 
Leriiis  were  included  in  the  radiative  dissipation  function.  In  contrast 
to  the  spherical  moilol,  tlie  effect  of  including  radiation  reaction  teri.is 
was  to  Increase  tlie  energy  loss. 
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